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Abstract 

We study the semiclassical dynamics of a spinor condensate with the magnetic dipole-dipole in- 
teraction included. The time evolution of the population imbalance and the relative phase among 
different spin components depends greatly on the relative strength of interactions as well as on the 
initial conditions. The interplay of spin exchange and dipole-dipole interaction makes it possible to 
manipulate the atomic population on different components, leading to the phenomena of sponta- 
neous magnetization and Macroscopic Quantum Self Trapping. Simple estimate demonstrates that 
these effects are accessible and controllable by modifying the geometry of the trapping potential. 

PACS numbers: 03.75.Kk, 03.75.Mn, 03.75.Nt, 05.30.Jp 
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I. INTRODUCTION 



One of the current experimental advances of considerable importance in the context of 
ultracold atoms is the demonstration of multicomponent Bose-Einstein condensate (BEG). In 
^articular several groups have created spinor condensates of ^^Na atoms [1] and ^'^Rb atoms 
3] by transferring spin polarized condensates into far-off-resonance optical dipole traps, 
where the spin degree of freedom becomes active and the spinor nature of the condensates is 
manifested. The spin-mixing dynamics and its dependence on the magnetic field have been 
investigated in detail for both spin-1 and spin-2 condensates in a number of theoretical yl 

riQ 

and experimental works Since the spin degree of freedom becomes accessible in an 

optical trap, the magnetic dipole-dipole interactions (MDDI) which arise from intrinsic or 
induced field of magnetic dipole moment 0, 0] should be taken into account in addition to 
the isotropic s-wave contact interaction. Due to their long-range and vectorial characters, 
the MDDI may largely enrich the variety of phenomena observed in condensates. Therefore 
this opens up new directions of research, such as the ground-state structure of spin-1 dipolar 



condensates in a single trap 



as well as the ground-state magnetic properties of spinor 



BEG confined in deep optical lattices , where only the dipolar interactions between 

condensates at different sites are considered. Recently, a beautiful experiment in Stuttgart 
[lol | has demonstrated BEG in a gas of chromium atoms ^^Cr, the MDDI of which is a factor 
of 36 higher than that for alkali atoms. This achievement makes possible the studies of the 
anisotropic long-range interactions in degenerate quantum gases. 

The present paper exploits semiclassical dynamics of spin-1 dipolar condensates in a sin- 
gle trap and shows the time evolutions of the population imbalance and the relative phase 
among different spin components 0, . This study may help us gain some insights into the 
properties of this dipolar spinor condensate, such as quantum phase diffusion, spontaneous 
magnetization jl^ and Macroscopic Quantum Self Trapping (MQST). Especially, MQST 
known as a novel nonlinear effect has already been pred icted theoretically and very 
recently observed experimentally in double- well trap [l4j and periodic optical lattice (l5l |. 
The self trapping effect sustains a self-maintained interwell population imbalance during the 
nonlinear tunneling process. Here we will show that for a dipolar condensate the interplay 
between spin-exchange and MDDI gives rise to the intercomponent MQST naturally. More- 
over, due to the adjustability of the MDDI, the experimental observation of the MQST effect 
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in dipolar spinor condensates can be expected. 

The paper is organized as follows. In Sec. II we describe the dipolar spinor BEC model 
and derive the equations of motion of semiclassical dynamics in terms of symmetry. Sec. Ill 
is devoted to the properties of the equations of motion such as spontaneous magnetization 
and spin-mixing dynamics. The time evolutions of the population imbalance and the relative 
phase among different spin components are investigated and quantum tunneling and self 
trapping among different spin components are examined carefully. Finally, a brief summary 
is given in Sec. IV. 



II. THE MODEL AND THE EQUATIONS OF MOTION 

We consider a spin F = 1 dipolar condensate with bosons. In second quantized form, 
the Hamiltonian Htot subject to both spin-exchanging collisions Hgp and MDDI Hdd reads 



Htot — Hsp + Hdd 



with 



Hsp = J dr^i (r) + Vt (r) I (r) 

+ y y dYi)l (r) (r) (r) i)f3 (r) 

+ dv^i (r) (r) F,,, ■ (r) (r) (1) 



and 



Hdd = ^ [ dr I dr'- 



2 J J |r — r'l 
X [^i (r) (r') F,„, ■ Ffs^,i^fs> (r) V^,, (r') 

- 3g (r) (r') (F„,, ■ e) (F^^, ■ e) i^^, (r) t^,, (r') ] , (2) 

where ijja (r) {a = 0, ±1) are the field annihilation operators for an atom in the hyperfine 
state \F = l,mp = a). The hyperfine spin F of one single atom is expressed in the spin-1 
matrix representation, and M is the mass of the atom. A summation over repeated indices 
is assumed in Eqs. ((H) and The external trapping potential Vr (r) is spin independent 
for a far off-resonant optical dipole trap which makes atomic spin degree of freedom com- 
pletely accessible. The two coefficients Cq = 47r^^ (oq + 202) /3M, C2 = inli^ (02 — Oq) /3M 



characterize the density- density and spin-spin colhsional interactions, respectively, with a/ 
(/ = 0, 2) being the s-wave scattering length for two spin-1 atoms in the combined symmet- 
ric channel of total spin /. And the dipolar interaction parameter is q = fiQu'^^gp/ATr with 
Qf being the Lande (7-factor, fis the Bohr magneton, e = (r — r') / |r — r'| is a unit vector. 
For the two experimentally created spinor condensates (^^Na and ^''Rb), we have \c2\ <^ Cq. 

For a condensate with its symmetry axis chosen to be along the quantization axis z (which 
happens to be the most experimentally relevant cases), Hm takes a very simple form as 
shown in Ref. ^j. Under the single mode approximation (SMA), namely, ipa (r) ~ 
where (r) is the spin-independent condensate spatial wave function, is the annihilation 
operator for mp = a component, the total Hamiltonian apart from a trivial term is given as 

a a,l3 

+ U2{o\a\aiai + a^_-Ji^_-Ji^ia^i — 2ajal]^aia_i 

-|- 2a\alaiao + 2aLiaoa_iao + 2a|)ajaia_i + 2a\a}_^aoao) 

+ Ud{2a\a\aiai + 2a^a^a_ia_i — 4a\a^_^aia_i 

— 2a|aJaiao — 2a^j^aJa_iao — 2aoaJaia_i 

— 2a\a^_^aQaQ + a\ai + a^_-fi_i — 2aJao)) (3) 

where s = J drip* (r) [— /2M + Vt (r)]0 (r) is assumed to be of the same value for all 
components, f/0,2 = (co,2/2) / |0 (r)|^(ir characterizes the spin-changing collisions and Ud = 
(crf/4) / / dvdr' |0 (r)|^ |0 (r')|^ (1-3 cos^ 6^) / |r - r'|^ denotes the MDDI with 6^ being the 
polar angle of (r — r'). Before we proceed to examine the semiclassical dynamics of Htot-, we 
would like to emphasize that H^p is invariant under spin rotation, i.e., it possesses a S'0(3) 
symmetry in spin space Q]. The presence of the MDDI breaks this symmetry into 5*0(2), 
which means Htot has only axial symmetry in spin space as we chosen the quantization axis 
along z. In this sense the MDDI plays a similar role to that of an external magnetic field. 
In the mean-field theory the condensate is usually considered to be in a coherent state 

1^) = exp I -i ^ ) exp [^Zadl ) |0), (4) 



a 



where |0) is the vacuum state. The complex numbers Za are nothing but the macroscopic 
wave functions for the atoms in the hyperfine level \F = l,mF = a) with population N^^ 



and phase 6a, i.e., 

Za = ^ae-''-. (5) 
The time-dependent variational principle of the system 

5S = 5 j ih{~z\'^) - {~z\H\'z)dt = Q (6) 

gives rise to the Hamiltonian equations of motion in complex coordinates 



a 



with 

Ho (iVa, Oa) = {e + U,)N- 3U,No + U^N^ + {U^ + 2Ud) {N, - N_,f 

+ 2 (f/2 - Ud) iVo[iVi + iV_i + 2^Njr[cos {260 -9,- e_,)]. (8) 

Note that the dynamics governed by Hamiltonian (j2)) conserves the total atom numbers N = 
X]q! and the total hyperfine spin of the condensate in the direction of the quantization 
axis / = (A^^i — A^_i) h. In terms of the following canonical variables 

^2 = 00- ^^i^ and = fiVo - i (iVi + iV-i) (9) 
the Hamiltonian (jHJ becomes cyclic in the coordinates and <^3 

(^^1, ^^2, ^^3, ^2) = + Uonj - 3f/d^]2 

+ 4 (f/2 + 2Ud) nl + 2{U2- Ud) 

X [Qr]i+fi2') {^^1-^2 



+ w Qr^i - vt^ - Ani (^VLi + ^2^ cos 2^2]- 



(10) 



Two important consequences follow as a result of the cyclic coordinates (pi and (p^. Firstly 
the mean value of the total number of atoms N = VLi and that of the total hyperfine spin 
of the condensate I = 2Ml^ are constants of motion. The dynamics, on the other hand, 
involves only one pair of variables {<y92, ^2} as 

. dHo ■ dHo 

^2 = T^, ^^2 



dil2 dip2 
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Define = Vti/N {i = 1, 2, 3) and rj = ^{^ - 6)^ - 4^1, we obtain 

= 4iV {U2 - Ud) r/ Q + 6^ sin2yP2 

+ 2N (f/2 - f/.) (i - (^^^ - ^) + (11) 

That means not only ^2 but also ^3 are involved in the dynamics. 

The variables ^2 and are canonically conjugate in a classical Hamiltonian 

+ 2 (f/2 - Ua) Nt] Q + cos 2(^2 - 3f/rf6- (12) 

In a simple mechanical analogy, describes a nonrigid pendulum. The Bose Josephson 
junction (BJJ) tunneling current between different spin components is given by 

^2N 



l-BJJ 



77 + ^2 sin 2(^2, (13) 



where Jq = 2N'^ (f/2 — Ud)- It differs from BJJ tunneling current of two weakly linked BEC 
in a double- well potential in its further nonlinearity in ^2 |il3i] . 



III. PROPERTIES OF THE EQUATIONS OF MOTION 
A. Spontaneous magnetization 

We firstly study the equilibrium configurations of the system. They are determined by 
the classical equations of motion after setting the time derivative terms to zero 

= r?Q + 6)sin2y;2 (14) 

= A ± [26 - \ 

^(tM£6^Meos2,2]. (15) 

We have defined a dimensionless parameter A = 2\U2-Ud\N characterizing the relative inter- 
action strength of the MDDI. In Eq. |T3|) plus sign corresponds to the case of {U2 — Ud) > 0, 



while minus sign corresponds to {U2 — Ud) < (In the following we take the former case as 
an example). From the definition of rj we know that the constant of motion ,^3 takes values 
in the interval — | < ^3 < | and the dynamic variable ^2 in the interval — | < ^2 < f ~ 2 |^3| . 
In the following discussion about the solutions of the equilibrium equations we separately 
consider three cases in order to see more clearly how these solutions depend on the relative 
phase ip2- 

(1) Equilibrium configuration with cos 2^92 = 1 (or 2ip2 = 2fc7r). 

The equilibrium value of ^2 is thus given by equation (fT3j) with cos 2(f2 = 1, which has only 
one solution in the interval — 00 < A < 1 + \J I — (2^3)^. When A ^ — 00 the equilibrium 
value of ^2 approaches its upper boundary | — 2 |^3|. In this case the fractions of atoms 
Ua = Na/N (a = 0, ±1) occupying three hyperfine states are ni = |^3| + ^3, no = 1 — 2 |^3| , 
= 1^3 1 — ^3 . On the other hand when A = 1 + ~ (2^3)^) ^2 is at the lower boundary 
— |, where the fractions of atoms occupying the hyperfine states are ni = ^ (1 + 2^3) , no = 0, 
n_i = I (1 — 2,^3) . Besides, when ^3 = the equilibrium value of ^2 is G = | (| ~ ^) in the 
interval —2 < A < 2 and the occupation fractions of three hyperfine states are respectively 
no = 1(1 - A) , m = n_i = i (1 + A) . 

(2) Equilibrium configuration with cos2(y92 = —1 (or 2ip2 = {2k + 1) tt). 

Now the equilibrium value of ^2 is given by equation (jl5|) with cos 2^92 = — 1- Similarly it 
has again only one solution in the interval A > 1 — ^1 — (2^3)^. When A = 1 — ^1 — (2^3)^, 
,^2 is at the lower boundary — |. In particular, when ^3 = the solution only exists for 
vanished dipole-dipole interaction Ud- 

(3) Equilibrium configuration with cos2(/32 = (or Equilibrium configuration which does 
not depend on the phase ip2)- 

In this case there exist two solutions which are independent of the value of A. One solution 
appears at the lower boundary ^2 = — | and the occupation fractions are rii = n_i = |, no = 
0. The other one appears at the upper boundary ^2 = f — 2 |^3| and the occupation fractions 
are rii = |,^3| + ,^3, no = 1 — 2 |,^3| , n_i = 1.^3 1 + ^3- Interestingly when .^3 = the occupation 
fractions become ni = n_i = and no = 1, corresponding to the "polar" state of the 
condensate. 

The mean field theory predicted a polar ground state of spinor BEG for ^^Na in which the 
atoms interact antiferromagnetically with each other in the absence of the MDDI Q]. The 
intrinsic magnetic moment of particles may, however, contribute together to form a Bose- 
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Einstein Ferromagnetism (BEF) Spinor bosons carry magnetic moments /i and M = 

filipif represents the magnetization. Spontaneous magnetization means that M remains 
finite even if the external magnetic field B vanishes, namely, the ground state of spinor 
BEC is always ferromagnetic state. Our calculations prefer the latter scenario. Specifically, 
above results suggest that for most given values of A and initial relative phase among the 
three components, when B = and U2 > 0, we have M 7^ , i.e., the system is magnetized 



121, one can 



spontaneously once the dipolar spinor BEC is realized. As mentioned in Ref. 
provide a direct way of confirming spontaneous magnetization experimentally. Considering 
that both the sign and the ma gnit ude of the MDDI coefficient Cd depend on the geometric 
shape of the condensate P. Il^ IitI. IfS^ . we can conveniently adjust the trap aspect ratio to 
manipulate the magnetic property of its ground state. 



B. Spin-Mixing Dynamics 

A key feature of dipolar spinor BEC is that besides the usual two-body repulsive hard- 
core interactions, there also exist spin-exchange interaction and MDDI which lead to spin 
mixing within the condensates. Population can be transferred from one spin state to another 
under internal nonlinear interactions without the presence of external fields. Insight into the 
complex dyiiamics of our system can be gained by employing the method of action-angle 
variables 1^1 . In the subsequent section we present some numerical results and show the 
time evolutions of the population imbalance and the relative phase among different spin 
components for two ) ^3 = and b) ,^3 7^ 0, i.e., whether total hyperfine spin of the 

condensate is zero or not. Time has been rescaled in units of 2 \U2 — Ud\ N/h in Figures 1 
and 2. 

Figure 1 shows solutions of the population imbalance and relative phase of Eqs. 
with initial conditions ^2(0) = 0.12, ^3(0) = 0.25, <f2{0) = f for the relative interaction 
parameters A = 2\U2-Ud\N ~ 0.00,0.40,0.60 and 0.75, respectively. The left column exhibits 
the time evolutions of the population imbalance between and ±1 components. Josephson 
sinusoidal oscillations are observed as A increases in Figs, la, lb. In Fig. Ic, there is a 
critical transition point for A = Ac = 0.60. Increasing A further for example to A = 0.75, the 
population of each equivalent component oscillates around a nonzero time averaged value, 
which gives a net population imbalance (^2 (t)) > 0. This is closely related to the MQST 
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phenomenon. Simultaneously, the right column shows the time evolutions of the relative 
phase between and ±1 components and indicates that as A increases (y92 (t) varies from a 
monotonically increasing function of time to a periodically oscillating function of time. In 
the nonrigid pendulum analogy, this corresponds that the motion of our system is turned 
from "running-phase modes" (Figure la) into "vr-phase modes" (Figures lb- Id) or from a 
rotation into a vibration. Here the definitions of running-phase modes and vr-phase modes 
are the same as those for two weakly coupled BEC [l^ . 

This critical behavior depends on = A^ [^2 (0) , ^3 (0) , ip2 (0) = f ] , as can be easily 
found from the energy conservation constraint and the boundness of the tunneling energy 
in Eq. (fT^ . In fact, the value (^2 (t)) < is inaccessible at any time if 

A>A.4-6(0) + ^v/^^^ 



3 ' 96 (0) 



1 + 36 (0))V(2-36(0))^-366M0) 



96 (0) 

When ip2 (0) = the critical parameter is 



(16) 



. _ 1 , ... Vl-9a(0) 
^^-3-^'^^^ 96 (0) 



1 + 36 (0))v'(2- 36(0))' -36^1(0) 



(17) 



96(0) 

and A < Ac marks the regime of MQST. Specifically for 6 (0) = 0.12, 6 (0) = 0.25, the 
critical parameter is Ac = —0.18. We take A = —0.20 and see MQST does set in. 

In Figure 2 we show solutions of Eqs. (fTT|) with initial conditions 6 (0) = 0.12, 6 (0) = 0, 
ip2 (0) = for parameters A = 0.75, 0.60, 0.43 and 0.00, respectively. The left column shows 
again the time evolutions of the population imbalance between and ±1 components and 
indicates that 6 (t) is always a periodic function of time as A decreases. For A = Ac = 0.43 
the population difference is self locked to the initial value, which serves as another sign of 
the MQST phenomenon. The right column shows the time evolution of the relative phase 
between and ±1 components and indicates that as A decreases {t) is always a periodic 
function of time around its mean value {ip2 (t)) = 0. The dynamics corresponds to "zero- 
phase modes". Moreover we observe that MQST occurs when 

A < Ac = ^ - 26 (0) , (18) 



while for (^92 (0) = vr MQST never happens. 

Therefore it is obvious that the spin-mixing dynamics depends on the relative interaction 
parameter A and is also sensitive to the initial occupations and relative phases of the three 
components, which can be adjusted by engineering Raman pulses js^. In practical experi- 
ments there are usually two different ways to achieve MQST. In Figures 1 and 2, ^2 (0) and 
(p2 (0) are kept constants while A varies (by changing the geometry of condensates). On the 
other hand, one can calibrate the initial values of the population imbalance ^2 (0) with a 
fixed trap geometry (i.e., A remains constant) and <^2 (0) 15 1. 

In order to further characterize the evolution of the system we summarize the full dynamic 
behavior of Eq. (jllj) in Figure 3 that shows the ^2 (^)-V^2 (^) phase portrait with constant 
energy contour. The distinction between the two dynamic regimes - nonlinear Josephson 
tunneling and MQST - becomes more apparent in Figure 3. In the regime of Josephson 
oscillation the dynamic variables follow a closed phase space trajectory, while in the self- 
trapping regime they follow an open trajectory with an unbounded phase. 

We are inspired by the great expectation that tuning the contact interaction between 
cold atoms close to zero by Feshbach resonance 0] will make the MDDI more prominent 
or even the dominant interaction [10]. Based on the experimental values of the s-wave 
scattering lengths for ^^Na, oq = (50.0 ± 1.6) as and 02 = (55.0 ± 1.7) as with as being the 
Bohr radius 23|, the ratio of coefficients q andjc2| can be shown to be 0.007, while for 
^*'^Rb(ao = 101 . Sor and 02 = 100. 4as) it is 0.1 p]. Evaluation with a simple variational 
wave function j2J| gives U2N ~ QnK for a sodium condensate up to = 5 x 10^ atonis 
confined in an optical dipole trap with a very sma 



1 trapping volume (10 cm ) [jj, 12 



On the other hand, as argued in Refs. ^,0, 10, 17, |l^, both the sign and the magnitude 
of the MDDI Ud can be greatly tuned with trapping geometry. As has been shown Q], 
Ud/U2 depends on a monotonically increasing function of the condensate aspect ratio k, 
bounded between —1 and 2. This clearly shows the possibility of adjusting the MDDI very 
close to spin exchange interaction, i.e. f/d/f/2 — 1. At this point our parameter A may 
take values in a large scale and the manifestation of MQST is within reach with current 
technologies. We point out although for all of the alkali-metal atomic condensations the 
MDDI is rather weak compared to the contact potential, the experimental achievement of 
BEG with transition-metal chromium ^^Cr provides us the hope because the MDDI here is 
36 times stronger than that of alkali atoms. For these reasons, a degenerate quantum gas 
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with adjustable long- and short-range interactions can be experimentally realized in near 
future. By altering the strengths of two kinds of interactions and the initial conditions, the 
nonlinear tunneling dynamics of dipolar spinor BEC consequently sustains a self-maintained 

population imbalance: a novel MQST effect. Such a scheme would avoid the difficulty of 
realizing experimentally the double-well magnetic trap. 

IV. SUMMARY 

In summary we have described a semiclassical treatment of the spin-1 dipolar spinor 
condensates. As a result of the conservation of atom numbers and total hyperfine spin of 
the condensate, the classical equations of motion are derived and discussed in a similar way 
as in double-well BJJ. It is demonstrated that spontaneous magnetization and spin-mixing 
dynamics depend on both the spin-exchange interaction U2 and the MDDI Ud through 
the ratio 2\U2-Ud\N - "^^^ initial population imbalance, the relative phase among the three 
components of the condensate as well as the total hyperfine spin of the system all play 
important roles in the semiclassical dynamics. Finally we have indicated the possibility 
of using dipolar spinor condensate as a platform for practical manipulation of MQST and 
Josephson oscillation. 
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Figure Captions 

Figure 1: Spin-Mixing Dynamics: the population imbalance ^2 and relative phase (^2 versus 
time for initial conditions 6 (0) = 0.12, ^3 (0) = 0.25, (0) = f and A = (a) , 0.4 (6) , 0.6 
(c) and 0.75 (d). Time is rescalcd in units of 2\U2 — Ud \ N/h. 

Figure 2: Same as in Figure 1 for initial conditions ^2 (0) = 0.12, ^3 (0) = 0, (p2 (0) = and 
A = 0.75 (a) , 0.6 (6) , 0.43 (c) and (d) . 

Figure 3: Constant energy contours in a phase-space plot of population imbalance ^2 versus 
phase difference (p2- The upper panel: ^3 (0) = 0.25 and A = 0.6 (c2 > 0). The lower panel: 
^3 (0) = and A = —0.2 (c2 < 0). Dark areas mean lower energy. 
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